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Abstract Nanoclusters of nickel sulphide (NiS)„ and iron 
sulphide (FeS)^ for n = 3-5 are studied using B3LYP 
exchange correlation function with 6-31G as a basis set. 
The structural stability of different isomers of NiS and FeS 
is analysed with the optimized energy, binding energy and 
vibrational studies. The electronic properties of isomers are 
discussed in terms of HOMO-LUMO gap, ionization 
potential, electron affinity, and embedding energy of dif¬ 
ferent clusters. Based on the calculated energy, planar ring, 
linear ladder and bipyramidal cube of NiS and FeS na¬ 
noclusters are found to be more stable. The increase in 
number of atoms in the clusters leads to increase in its 
stability. The dipole moment is high for planar rhombus 
and linear ladder structures of NiS and FeS nanoclusters. 
The ionization potential and electron affinity are high for 
planar ring structure of NiS and FeS clusters. The value of 
energy gap of linear ladder NiS nanocluster and cube and 
bipyramidal cube structures of FeS is found to be low. The 
binding energy for cube structure of NiS and FeS clusters is 
found to be high. Hexagonal ring structure of NiS and FeS 
clusters has low embedding energy value. The other 
parameters such as dipole moment and point symmetry are 
also discussed. The structural stability and electronic 
properties will provide an insight for experimentalist to 
tailor new materials that have potential importance in 
engineering applications. 
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Background 

Chalcogenides are an important class of materials which 
have significant electrical properties, optical properties and 
chemical characteristics. Chalcogen includes sulphur, 
selenium and tellurium elements. The sulphur compounds 
have a wide range of properties which attract the scientific 
community, predominantly in thin film technology and 
nanoparticle synthesis. The applications of chalcogenide 
materials include a variety of chalcogenide glasses, infra¬ 
red sensors, solar energy conversion and window layer 
[1-6]. Various methods are utilized for the synthesis of 
binary, ternary and quaternary chalcogenides. 

In the family of chalcogenide nanomaterials, nickel 
sulphide (NiS) is one of the important materials which 
find its potential application in photoconduction [7]. Han 
et al. [8, 9] reported NiS as an attractive cathode 
material for rechargeable lithium batteries and charge- 
discharge mechanism of NiS as a cathode material. 
Fernandez et al. [10] studied ZnS-NiS-CuS optical filters 
with wide-range solar control characteristics and as a 
catalyst [11]. Different morphologies of NiS can be 
prepared by various methods: Shen et al. [12] reported 
phase-controlled synthesis of nickel sulphide nanorods; 
Wang et al. [13] synthesized NiS nanotubes using self- 
assembly route in AAO templates; Chen et al. [14] 
prepared NiS nanotubes and nanoneedles derived from 
rolled nanosheets from microemulsion; and Huang et al. 
[15] reported NiS with a core-shell structure. NiS can be 
synthesized using solvothermal method [16] and hydro- 
thermal process [17]. 
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Fig. 1 Structures of NiS 
nanoclusters 



(a) NiS - Hexagonal ring 



(b) NiS - Planar rhombus 


1.9 A 



(c) NiS - cube 


(d) NiS - Planar ring 



(e) NiS - linear ladder 


(f) NiS - bi pyramidal cube 


Table 1 Energy and dipole moment of geometrically optimized NiS nanoclusters 


Size 

Model 

Energy (Hartrees) 

Dipole moment (Debye) 

Symmetry 

3 

Hexagonal ring 

-5,718.76 

1.1189 

C s 

3 

Planar rhombus 

-5,718.44 

3.2665 

C s 

4 

Cube 

-7,623.87 

0.3002 

Cl 

5 

Planar ring 

-9,531.34 

0.2135 

Cs 

5 

Linear ladder 

-9,531.09 

3.7378 

C s 

5 

Bipyramidal cube 

-9,530.06 

1.4671 

Cl 
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Table 2 Energy gap and density of states of NiS nanoclusters 

S. No. Cluster type Energy gap(eV) NiS 


1 


Hexagonal ring 1.04 



Energy (eV) 


2 Planar rhombus 1.79 


3 


Cube 


1.5 
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Table 2 

continued 



S. No. 

Cluster type 

Energy gap(eV) 

NiS 


4 


Planar ring 


1.22 


5 Linear ladder 0.82 


6 


Bipyramidal cube 1.53 



Energy (eV) 


w 


Springer 






















































































































































































J Nanostruct Chem (2014) 4:87 


Page 5 of 16 87 


Iron sulphide (FeS) has importance in the function of 
electron-transfer groups in biological applications such as 
enzyme catalysis, protein folding and metabolic regulation, 
photosynthesis and respiratory chain. Esrafili et al. [18], 
studied electronic structure of small Fe m S„ nanoclusters; 
Henderson et al. [19] reported FeS-based materials for 
groundwater remediation; Ritchie et al. [20] lithium-ion/ 
iron sulphide rechargeable batteries; and Gomes et al. [21] 
synthesized FeS by electrodeposition method. 

The motivation of the present work is to tailor proper 
structure of NiS and FeS with suitable method that will 
enhance its structural and electronic properties with 
improved electrical, optical and chemical properties. With 
this as motivation, the survey was conducted using 
Thomson Reuters “web of science” database, only less 
work has been reported in NiS and FeS. The reported work 
mainly concentrates on synthesis and characterization of 
nickel sulphide and iron sulphide. The quantum chemical 
studies on these chalcogenide materials will fine-tune the 
structural and electronic properties which will be suitable 
for its engineering application. Using the International 
Centre for Diffraction Data (ICDD) data as reference, the 
nanoclusters are constructed; NiS is taken from card 
number: 89-7141 and FeS is taken from card number: 
89-9021. The density functional theory (DFT) is a better 
approach to study the structural and electronic properties 
for chalcogenide materials. The electron movement is 
confined in nanoclusters which resemble zero-dimensional 
material. The structural and electronic property mainly 
depends on the geometry of the clusters; DFT method 
effectively explores the electronic properties of nano¬ 
structures. With this as the objective, some of the possible 
six realistic structures for each group are studied using the 
NWChem package and the results are analysed and 
reported. 

Methods 

The realistic clusters of NiS and FeS of different structures 
were completely optimized using NWChem package [31]. 
The optimization in this work is based on B3LYP [32-35] 
exchange correlation which completely exploits sulphide 
clusters with 6-31G as a basis set. Since, the atomic 
number of nickel, iron and sulphur are 28, 26 and 16, 
respectively; 6-31G is the suitable basis set which will 
optimize these clusters [36-38] and this is the reason 
behind the selection of this basis set. In the quantum 
chemical methods, the pseudo potential approximation will 
replace the complex field of bound electrons in the atoms 
and the effective potential of the nucleus is a modified 
potential term. NWChem utilizes Fock matrix elements, 
which in the case of Hartree-Fock calculation with the 



Cluster 

Fig. 2 Variation of IP and EA with cluster type for NiS clusters 

local basis set called Fock build is adopted; a similar 
algorithm of Kohn-Sham matrix elements is used to 
evaluate the computation in DFT calculations [39]. The 
convergence study was carried out with 10~ 5 eV with 
respect to energy. 

Results and discussion 

The realistic possible structure of sulphide clusters of NiS 
and FeS is designed in six different types. The hexagonal 
structure resembles bees hive like structure. The planar 
rhombus structure has rhombus structure which is sepa¬ 
rated by square. The cube structure has a cube-like 
arrangement. Planar ring structure has five transition metal 
atoms and five sulphur atoms which form a closed ring-like 
loop. In the linear ladder arrangement, all the atoms form a 
ladder-like structure. In the bipyramidal structure, the cube 
is projected with the pyramid at two faces. These are the 
different types of nanoclusters which are studied in the 
present work. The energy, dipole moment, point symmetry, 
HOMO-LUMO gap, ionization potential (IP), electron 
affinity (EA), binding energy (BE), embedding energy (EE) 
and vibration analysis of different sulphide clusters are 
analysed and reported. 

Structures of (NiS)„ 

The possible geometrically optimized structures of (NiS)„ 
are shown in the Fig. 1. The hexagonal ring and planar 
rhombus have the energy of —5,718.76 and —5,718.44 
Hartrees, respectively, with the dipole moment of 1.118 
and 3.266 Debye, respectively. Due to closed structure of 
hexagonal ring the dipole moment is comparatively lower 
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Fig. 3 Structures of FeS 
nanoclusters 


1.9 A 




(c) FeS - cube 

1.9 A 



(e) FeS - linear ladder 


(b) FeS - Planar rhombus 




(d) FeS - Planar ring 


(f) FeS - bi pyramidal cube 


Table 3 Energy and dipole moment of geometrically optimized FeS nanoclusters 


Size 

Model 

Energy (Hartrees) 

DM (Debye) 

Symmetry 

3 

Hexagonal ring 

-4,984.93 

1.5816 

C s 

3 

Planar rhombus 

-4,984.62 

4.1655 

C s 

4 

Cube 

-6,645.53 

0.8781 

Ci 

5 

Planar ring 

-8,308.3 

0.9713 

C s 

5 

Linear ladder 

-8,308.02 

5.1115 

C s 

5 

Bipyramidal cube 

-8,307.09 

1.4485 

Cl 


than planar rhombus. The cube structure has the energy of 
—7,623.87 Hartrees, cube structure has low value of 
0.300 Debye which is due to the ordered structure. The 
planar ring, linear ladder and bipyramidal cube have the 


energy in the order of —9,530 Hartrees. In these structures 
with (NiS) 5 , the dipole moment is large for linear ladder 
structure whereas, the other two are almost symmetrical 
and have the low value of dipole moment. The obtained 


w 


Springer 















J Nanostruct Chem (2014) 4:87 


Page 7 of 16 87 


Table 4 Energy gap and density of states of FeS nanoclusters 

S. No. Cluster type Energy gap(eV) FeS 


1 


Hexagonal ring 1.66 


2 Planar rhombus 2 


3 


Cube 


1.35 




— DOS spectrum 

- Occupied orbitals 

- Virtual orbitals 




—1 — LI -III II II.- II II II - l»l - 1— LI - JJ - 11 H I 

-10 -5 0 5 10 

Energy (eV) 



Energy (eV) 


energy and dipole moment of geometrically optimized 
clusters are tabulated in Table 1. 

The HOMO-LUMO gap gives the electronic properties 
of NiS nanoclusters [22, 23]. Usually a low value of gap 


implies that the electron can easily move from the occupied 
level to the virtual level. Among all the different structures 
of NiS clusters, a low value is seen for linear ladder which 
is more suitable for electronic applications. In contrast, a 
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Table 4 continued 



S. No. Cluster type 

Energy gap(eV) 

FeS 


4 


Planar ring 


1.62 



Energy (eV) 


5 Linear ladder 1.61 


6 


Bipyramidal cube 1.33 



high value of gap is observed for planar rhombus. In the 
planar rhombus Ni and S atoms are arranged symmetrically 
that gives rise to the high value of HOMO-LUMO gap due 
to overlapping of unfilled shells in d and p orbitals of Ni 
and S atoms. The values of HOMO-LUMO gap of NiS 
nanoclusters are tabulated in the Table 2. 


The high value of IP implies that it is difficult to remove 
the electron from the cluster which will be more chemi¬ 
cally inert [24]. Among all the structures planar ring has 
the high value of IP due to the closed arrangement of atoms 
in the ring. EA refers the amount of energy change when 
the electron is added to the cluster [25, 26]; a low value 
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Cluster 


Fig. 4 Variation of IP and EA with cluster type for FeS clusters 

indicates less chemically reactive nature of the cluster. 
Low value of EA seen for cube structure is 2.63 eV. The IP 
and EA values of different clusters are shown in Fig. 2. 

Structures of (FeS)„ 

The structures of (FeS)„ are shown in the Fig. 3. The cal¬ 
culated energies and dipole moment of various FeS na¬ 
noclusters are shown in Table 3. The hexagonal ring and 
planar rhombus structure has energy of —4,984 Hartrees 
which has three Fe atoms and three sulphur in its clusters. 
In these structures a moderate dipole moment of 
1.581 Debye is observed for hexagonal ring, whereas a 


high dipole moment of 4.165 Debye is seen for planar 
rhombus. In the cases of planar ring, linear ladder and 
bipyramidal cubes, all have the energy in and around 
—8,308 Hartrees with different dipole moments due to the 
positions of atoms: planar ring has a low dipole moment of 
0.971 Debye, bipyramidal cube structure has a moderate 
value of 1.448 Debye and linear ladder has a high value of 
5.111 Debye. 

The high value of HOMO-LUMO gap is noticed for 
planar rhombus; the electron requires more energy to move 
from HOMO to LUMO level. The high value of gap in 
planar rhombus structures arise due to metal-metal bond¬ 
ing of Fe atoms which has unoccupied d shells. In the case 
of FeS clusters, all the clusters have a low value of gap 
which seems to have more metallic-like property. The 
HOMO-LUMO gap of FeS clusters are tabulated in the 
Table 4. 

The high value of IP around 5.58 eV is noticed for 
planar ring due to the closed loop. Comparing with the IP 
value of NiS structures, the FeS clusters have low values of 
IP. A low value of EA is observed for cube and bipyra¬ 
midal cube structures of FeS. The IP and EA variations for 
different clusters are shown in the Fig. 4. 

Binding energies and embedding energies of (NiS)„ 
and (FeS)„ clusters 

Binding energy is also one of the important parameter to 
analyse the structural stability of the nanoclusters [27]. The 
BE per atom of the nanoclusters can be calculated from the 
following Eq. (1). 


Table 5 Calculated binding energies of NiS nanoclusters Table 7 Calculated embedding energies of NiS nanoclusters 


Size 

Model 

BE (eV) 

Size 

Model 

EE (eV) 

3 

Hexagonal ring 

0.94 

3 

Hexagonal ring 

5.65 

3 

Planar rhombus 

2.39 

3 

Planar rhombus 

14.36 

4 

Cube 

4.84 

4 

Cube 

38.77 

5 

Planar ring 

0.76 

5 

Planar ring 

7.67 

5 

Linear ladder 

1.45 

5 

Linear ladder 

14.50 

5 

Bipyramidal cube 

4.25 

5 

Bipyramidal cube 

42.50 

Table 6 

Calculated binding energies of FeS nanoclusters 


Table 8 

Calculated embedding energies of FeS nanoclusters 

Size 

Model 

BE (eV) 

Size 

Model 

EE (eV) 

3 

Hexagonal ring 

1.59 

3 

Hexagonal ring 

9.55 

3 

Planar rhombus 

2.99 

3 

Planar rhombus 

17.98 

4 

Cube 

5.13 

4 

Cube 

41.11 

5 

Planar ring 

1.37 

5 

Planar ring 

13.76 

5 

Linear ladder 

2.14 

5 

Linear ladder 

21.49 

5 

Bipyramidal cube 

4.66 

5 

Bipyramidal cube 

46.63 
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Table 9 Vibrational spectrums of NiS nanoclusters 


S. No. Cluster type 

NiS 


Hexagonal ring 


Planar rhombus 


Cube 


Infrared Spectrum 



200 


400 


600 


1000 


Frequency 

Infrared Spectrum 



Infrared Spectrum 



BE = (n x E[TM] +n x E[S]-n x E[TMS])/n (1) 

where £{TM] is the energy of the transition metal, E[ S] is 
the energy of sulphur, £[TMS] is the energy of transition 
metal sulphide and n is the number of transition metal and 
sulphur atom in the cluster. The energy of individual atoms 
(TM and S) is calculated separately in the same basis set 
and the resultant energy values are taken into account to 
calculate £[TM] and E[ S]. BE plays a vital role in deter¬ 
mining the stability of the nanoclusters. 


Among all the clusters of NiS, BE for cube and bipy- 
ramidal cube was high which can be said to be more stable 
than that of other clusters; low BE is observed for hexag¬ 
onal ring and are tabulated in Table 5. Interestingly, in the 
case of FeS clusters also the cube and bipyramidal cube 
have the high BE. The BE values of FeS nanoclusters are 
tabulated in Table 6. The high value of BE for cube and 
bipyramidal cube ascends due to the similar kind of geo¬ 
metrical symmetry. Comparing the BE of all the clusters of 
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Table 9 continued 


S. No. 


Cluster type 


NiS 


Planar ring 


Infrared Spectrum 


Linear ladder 




Frequency 


6 Bipyramidal cube Infrared Spectrum 



Frequency 


NiS and FeS, high value is seen for FeS clusters. In NiS 
and FeS clusters, cube and bipyramidal cube have the high 
value of BE. The reason may be due to the perfect sym¬ 
metry of structure between the cube and bipyramidal cube. 

The other parameter which can be explored is the EE; it 
refers to the gain in energy to embed a foreign atom during 
the growth process in the cluster [28]. The EE of the cluster 
is given by relation as given in the Eq. 2, 


EE = £(TM)„+E(S) n -E [(TMS)J (2) 

where TM represents transition metals such as Ni, Fe, S for 
sulphur and n represents the number of transition metal 
atoms and sulphur atom present in the clusters. For the 
calculation of EE, the individual energy of TM and S atoms 
is calculated in the same basis set and is incorporated to 
find £(TM) and E( S). 
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The variation in the EE with the cluster size of NiS is 
shown in the Table 7. For NiS clusters, the bipyramidal 
cube structure of NiS cluster has a high value of EE which 
is difficult to add a foreign atom in this cluster. In FeS 
cluster, the EE is maximum for bipyramidal cube structure 
of FeS. The EE of FeS is tabulated in the Table 8. Among 
all the clusters, the EE of bipyramidal cube has the high 
value. 

Vibrational analysis of (NiS)„ and (FeS)^ clusters 
for n = 3-5 

The vibrational analysis provides the insight for the sta¬ 
bility of the nanoclusters [29, 30]. A particular cluster with 
a real positive frequency is said to be more stable. More IR 
intensity in the positive frequency leads to an increase in 
cluster stability. Tables 9 and 10 show the vibrational 
spectrum and mode assignment of NiS nanoclusters. Hex¬ 
agonal ring has the significant vibrational frequency at 
491.34 and 688.54 cm -1 with the intensity of 59.75 and 
27.36 km/mole both assigned to molecular stretching. In 
the case of planar rhombus, high IR intensity of 31.29 and 
89.59 km/mole is observed at the vibrational frequency of 
346.90 and 1,108.89 cm -1 assigned to molecular stretching 
and S-S stretching, respectively. Looking at the vibrational 
frequency of cube structure, the maximum intensity of 7.35 
and 21.12 km/mole for this structure is observed at the 
frequency of 951.55 and 1,047.37 cm 1 respectively mode 
assigned to molecular stretching. The prominent IR inten¬ 
sity of 26.18 and 42.05 km/mole is observed for planar ring 
at the frequency of 28.61 and 738.04 cm -1 respectively 
assigned to molecular stretching and S-Ni-S stretching. 


Linear ladder has the intensity of 12.09, 15.36 and 
20.27 km/mole at the frequency of 636.19, 655.15 and 
991.08 cm -1 which is mode assigned to molecular twist, 
molecular stretching and Ni-S-Ni stretching, respectively. 
Bipyramidal cube has the intensity of 20.16 and 21.32 km/ 
mole for the frequency of 634.23 and 1,070.82 cm -1 
assigned to molecular stretching and Ni-S-Ni stretching, 
respectively. 

The vibrational analysis of FeS is shown in Tables 11 
and 12. In the case of hexagonal ring, prominent IR 
intensity of 75.95 and 46.15 km/mole is observed at the 
frequency of 443.12 and 504.46 cm -1 for Fe-S-Fe 
stretching and Fe-S-Fe bending mode respectively. The 
planar rhombus structure has the vibrational frequency at 
741.44 and 1,118.84 cm -1 that has the IR intensity of 
12.57 and 61.50 km/mole for molecular stretching and S-S 
stretching respectively. The cube structure has the IR 
intensity of 17.24, 11.39 and 25.64 km/mole for frequency 
of 985.88, 1,110.02 and 1,137.97 cm -1 for the mode of 
molecular bending, molecular twisting and molecular 
twisting, respectively. Planar ring structure has the IR 
intensity of 39.33, 74.09 and 104.59 km/mole for the fre¬ 
quency of 570.5, 795.6 and 906.13 cm -1 assigned to 
molecular stretching, S-Fe-S stretching and Fe-S-Fe 
stretching, respectively. In the case of linear ladder, the 
vibrational frequency occurs at 671.1, 714.71 and 
752.79 cm -1 with IR intensity of 38.13, 330.66 and 
42.85 km/mole which is assigned to molecular stretching, 
molecular bending and Fe-S stretching, respectively. The 
bipyramidal cube has the vibrational frequency of 440.66, 
656.47 and 1,160.71 cm -1 with IR intensity of molecular 
twist, molecular stretching and Fe-S stretching. 


Table 10 Vibrational mode assignment of NiS nanoclusters 


S. No. 

Cluster size 

Model 

Frequency (cm x ) 

Intensity (km/mole) 

Mode assignment 

1 

3 

Hexagonal ring 

437.89 

23.63 

S-Ni-S stretch 




491.34 

59.75 

Molecular stretch 




688.54 

27.36 

Molecular stretch 

2 

3 

Planar rhombus 

346.90 

31.29 

Molecular stretch 




1,108.89 

89.59 

S-S stretch 

3 

4 

Cube 

951.55 

7.35 

Molecular stretch 




1,047.37 

21.12 

Molecular stretch 

4 

5 

Planar ring 

28.61 

26.18 

Molecular stretch 




738.04 

42.05 

S-Ni-S stretch 




1,049.28 

11.66 

Ni-S stretch 

5 

5 

Linear ladder 

636.19 

12.09 

Molecular twist 




655.15 

15.36 

Molecular stretch 




991.08 

20.27 

Ni-S-Ni stretch 

6 

5 

Bipyramidal cube 

634.23 

20.16 

Molecular stretch 




1,070.82 

21.32 

Ni-S-Ni stretch 
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Table 11 Vibrational spectrums of FeS nanoclusters 


S. No. Cluster type 

FeS 


1 Hexagonal ring Infrared Spectrum 



-0 200 400 600 800 


Frequency 

2 Planar rhombus Infrared Spectrum 



-0 200 400 600 800 1000 1200 

Frequency 

3 Cube Infrared Spectrum 



4 Planar ring Infrared Spectrum 
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Table 11 continued 


S. No. 


Cluster type 


FeS 


Linear ladder 


Infrared Spectrum 


Bipyramidal cube 




Table 12 Vibrational mode assignment of FeS nanoclusters 


S. No. Cluster size 

Model 

Frequency (cm x ) 

Intensity (km/mole) 

Mode assignment 

1 3 

Hexagonal ring 

443.12 

75.95 

Fe-S-Fe stretch 



504.46 

46.15 

Fe-S-Fe bend 



737.89 

15.62 

Fe-S bend 

2 3 

Planar rhombus 

741.44 

12.57 

Molecular stretch 



1,118.84 

61.50 

S-S stretch 

3 4 

Cube 

985.88 

17.24 

Molecular bend 



1,110.02 

11.39 

Molecular twist 



1,137.97 

25.64 

Molecular twist 

4 5 

Planar ring 

570.50 

39.33 

Molecular stretch 



795.60 

74.09 

S-Fe-S stretch 



906.13 

104.59 

Fe-S-Fe stretch 

5 5 

Linear ladder 

671.1 

38.13 

Molecular stretch 



714.71 

330.66 

Molecular bend 



752.79 

42.85 

Fe-S stretch 

6 5 

Bipyramidal cube 

440.66 

17.48 

Molecular twist 



656.47 

17.58 

Molecular stretch 



1,160.71 

34.72 

Fe-S stretch 

^ Springer 
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Conclusions 

Nanoclusters of (NiS) n and (FeS)„ for n = 3-5 are com¬ 
pletely optimized with B3LYP/6-31G as a basis set. The 
stability of the cluster increases when the cluster size 
increases. The dipole moment and point symmetry for all 
different types of isomers depend on the positioning of atoms 
in the cluster. The electronic properties are studied in terms 
of HOMO-LUMO gap, IP and EA of different isomers, 
among these clusters chemically reactive and inert clusters 
are identified and reported. The BE of different clusters are 
studied and reported. The embedding energy of transition 
metal sulphide clusters provides the insight for the inclusion 
of the foreign atom in the cluster. Vibrational analysis also 
reveals the stability of the clusters, the particular cluster with 
no imaginary frequency can be concluded as more stable. 
The structural and electronic properties of different isomers 
provide insight for tailoring new materials which find its 
potential importance in functional nanomaterials. 
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